Liu et al. BMC Plant Biology 201 1 , 1 1 :1 32 
http://www.biomedcentral.eom/1471-2229/1 1/132 



Plant Biology 



RESEARCH ARTICLE Open Access 



Mapping QTLs for oil traits and eQTLs for oleosin 
genes in jatropha 

Peng Liu, Chun Ming Wang", Lei Li, Fei Sun, Peng Liu and Gen Hua Yue"" 
Abstract 

Background: The major fatty acids in seed oil of jatropha, a biofuel crop, are palnnitic acid (CI 6:0), stearic acid 
(C18:0), oleic acid (C18:1) and linoleic acid (C18:2). High oleic acid and total oil content are desirable for jatropha 
breeding. Until now, little was known about the genetic bases of these oil traits in jatropha. In this study, 
quantitative trait locus (QTL) and expression QTL analyses were applied to identify genetic factors that are relevant 
to seed oil traits in jatropha. 

Results: Composite interval mapping identified 18 QTL underlying the oil traits. A highly significant QTL qCl 8:1-1 
was detected at one end of linkage group (LG) 1 with logarithm of the odd (LOD) 18.4 and percentage of variance 
explained (PVE) 36.0%. Interestingly, the QTL qC18:l-l overlapped with qC18:2-l, controlling oleic acid and linoleic 
acid compositions. Among the significant QTL controlling total oil content, qOilC-4 was mapped on LG4 a relatively 
high significant level with LOD 5.0 and PVE 11.1%. Meanwhile, oleosins are the major composition in oil body 
affecting oil traits; we therefore developed SNP markers in three oleosin genes Olel, Olell and Olelll, which were 
mapped onto the linkage map. Olel and Olelll were mapped on LG5, closing to QTLs controlling oleic acid and 
stearic acid. We further determined the expressions of Olel, Olell and Olelll in mature seeds from the QTL mapping 
population, and detected expression QTLs (eQTLs) of the three genes on LGs 5, 6 and 8 respectively. The eQTL of 
Olelll, qOlelll-5, was detected on LGS with PVE 11.7% and overlapped with QTLs controlling stearic acid and oleic 
acid, implying a cis- or trans-element for the Olelll affecting fatty acid compositions. 

Conclusion: We identified 18 QTLs underlying the oil traits and 3 eQTLs of the oleosin acid genes. The QTLs and 
eQTLs, especially qCl 8:1-1, qOilC-4 and qOlelll-5 with contribution rates (R^) higher than 10%, controlling oleic acid, 
total oil content and oleosin gene expression respectively, will provide indispensable data for initiating molecular 
breeding to improve seed oil traits in jatropha, the key crop for biodiesel production. 



Background 

Jatropha curcas is becoming one of the world's key 
crops for biodiesel production [1]. Oil containing a high 
amount of unsaturated fatty acid can find an application 
as biodiesel feed stock. To make the production of jatro- 
pha profitable and sustainable, genetic improvement of 
oil yield and quality is demanded. However, oil traits 
cannot be evaluated until the seeds are harvested and 
analyzed in laboratory, and detailed selective breeding 
has not been carried out. Meanwhile molecular breeding 
in jatropha has not been started due to lack of molecu- 
lar bases of economically important traits such as seed 
yield, seed oil traits, biotic or abiotic stress resistance. 
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Most economically important traits are quantitative 
and determined by many genes and gene complex 
where are described as quantitative trait loci (QTLs). 
Traditional methods of genetic improvement of quanti- 
tative traits have relied mainly on phenotype and pedi- 
gree information [2], which are easily influenced by 
environmental factors. To conduct marker assisted 
selection (MAS) for genetic improvement of oil yield 
and quality in jatropha, the molecular bases of seed oil 
traits need to be understood by identifying genomic 
regions that contain favorite loci, i.e. QTL analysis. QTL 
analysis has been performed to detect the genetic bases 
of important agronomic or physiological traits, providing 
valuable information for trait improvement. Genetic 
markers have made it possible to detect QTLs that are 
significantly associated with traits, and made selection 
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more effective [3]. Genetic response can be improved by 
including the QTLs in marker-assisted selection, which 
is a method of selection that makes use of phenotypic, 
genotypic and pedigree data [4]. Moreover, MAS for oil 
traits improvement will be much advantageous com- 
pared to traditional breeding because seed oil traits can- 
not be measured at early stage or in field. The use of 
DNA markers for selection in jatropha can greatly 
reduce breeding scale. By using MAS, decisions can be 
made at the nursery stage, regarding which individuals 
should be retained as breeding stock, and which should 
be removed. 

To conduct QTL analysis, most appropriate crosses 
need to be selected to generate sufficient genetic varia- 
tions both on DNA and phenotype levels. QTL analyses 
of total oil content have been made in a number of 
crops, including oilseed rape[5], soybean[6], maize[7], 
and sunflower [8]. Recent surveys have shown large var- 
iations in content and fatty acid composition of seed oil 
of Arabidopsis, suggesting populations derived from 
selected crosses will be useful for investigating these 
traits [9]. 

Diversity in gene expression is one of the mechanisms 
underlying phenotypic diversity among individuals and 
regarded as one of quantitative traits [10]. Analysis of 
determinants of candidate gene expression not only 
helps in understanding the mechanisms for phenotypic 
variation but also provides an approach to improve phe- 
notypes via the modulation of gene expression [10]. With 
advances in gene expression profiling, an approach 
named "genetical genomics" has been put forward to 
identify the determinants of gene expression [11]. This 
approach treats mRNA expression levels as quantitative 
traits in a segregating population and maps expression 
QTL (eQTL) that control expression levels in vivo. For 
almost any gene analyzed in a segregating population, 
eQTL analysis can identify the genomic regions influen- 
cing its expression level. eQTL that map to the same 
genetic location as the gene whose transcript is being 
measured generally indicate the presence of a cis-acting 
regulatory polymorphism in the gene (cis-eQTL). eQTL 
that map distant to the location of the gene being 
assayed most likely identify the location of trans-acting 
regulators (trans-eQTL) that may control the expression 
of a number of genes elsewhere in the genome. The 
genetical genomics approach has been employed for 
identifying eQTL regulating gene expression [10,12]. 

Recently, we established a first generation genetic link- 
age map of jatropha using 506 microsatellite and SNP 
markers covering 11 linkage groups [13], thus providing 
a necessary tool for a whole genome scan for QTLs and 
eQTLs affecting economically important traits including 
seed oil traits. Among the fatty acid present in the jatro- 
pha seed oil, linoleic acid (18:2), oleic acid (18:1), 



palmitic acid (16:0) and stearic acid (18:0) are dominant 
compositions. Oleic and linoleic acids are the major 
constituents of jatropha oil [14]. The breeding goal for 
jatropha seed oil trait improvement is to increase total 
oil content and oleic acid, and decrease palmitic content 
[15]. In this paper, we describe the genetic bases of 
these seed fatty acid composition and content traits 
through QTL mapping with a backcrossing population 
consisting 286 individuals. On the other hand, seed oil 
is stored in subcellular organelles called oil bodies. Pro- 
teome composition of the jatropha oil bodies revealed 
oleosins as the major component affecting oil traits [16]. 
Three jatropha oleosin genes, namely Olel, Olell and 
Olelll, were isolated [17]. Here, we developed SNPs of 
the three oleosin genes in the QTL mapping population, 
which were subsequently mapped onto the linkage map. 
We determined expression variations of the three genes 
in the QTL mapping population, conducted an eQTL 
analysis on oleosin gene expressions and provided new 
information for possible modulation of oleosin genes to 
improve oil traits in jatropha. 

Results 

Trait analysis 

Fatty acid composition, total oil content of jatropha 
seeds and gene expression levels of oleosin genes were 
measured in the QTL mapping population. The fre- 
quency distributions of the traits showed a continuous 
distribution (data not shown), revealing complex genetic 
bases of these traits. As expected for an interspecific 
cross, distribution of phenotypic values in the progeny 
showed bi-directional transgressive segregations for all 
traits (Table 1). C18:l in /. curcas is higher than in /. 
integerrima, while total oil content in /. integerrima is 
51.04%, much higher than in /. curcas. The data implied 
that /. integerrima germplasm could be applied for 
hybrid breeding to improve agronomic traits such as 
total oil content. 

Correlation analysis among these traits was per- 
formed (Table 2). CI 8:2 showed a significantly negative 
correlation with C18:l and C16:0. Especially the C18:l 
correlated with C18:2 with a high coefficient -0.962, 
implying that there could be common genetic factors 
affecting these two compositions. The expression levels 
of Olely Olell and Olelll showed a highly positive cor- 
relation with each other. Olel expression level was sig- 
nificantly correlated with C16:0. The correlation 
coefficients between expression levels of oleosin genes 
and total oil content were low but significant. The sig- 
nificant but low values of correlation coefficients 
implied genetic bases of fatty acid composition and 
total oil content were complex, and oleosin genes 
could be involved the multiple genetic factors affecting 
these oil traits. 
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Table 1 Descriptive statistics on phenotype data of the QTL mapping population and parents (J.curcas PZM16 and J. 
integerrima S001) 



Traits 


Mean 


SD 


Min 


Max 


PZM16(Mean ± SD) 


SOOKMean ± SD) 


CI 6:0 (%) 


12.6 


2.54 


7.95 


21.51 


18.88 ± 5.7 


7.89 ± 0.25 


CI 8:0 (%) 


6.28 


1.44 


2.89 


10.42 


6.02 ± 1.01 


5.13 ± 0.31 


CI 8:1 (%) 


37.72 


9.41 


18.44 


61.77 


42.42 ± 0.54 


30.83 ± 3.64 


CI 8:2 (%) 


43.39 


9.96 


20.57 


66.22 


32.7 ± 5.23 


56.14 ± 3.58 


Total oil content (%) 


35.4 


8 


13.3 


57.1 


30.59 ± 0.70 


51.04 ± 2.39 


Olel expression (AACj) 


-0.32 


3.63 


-9.89 


8.97 


0 


-0.42 


Olell expression (AACj) 


2.13 


3.78 


-6.03 


11.33 


0 


2.54 


Olelll expression (AACj) 


-3.4 


4.05 


-12.1 


3.63 


0 


0.98 



QTL and eQTL mapping 

The linkage map, covering 663.0 cM of the genome, 
converged into 11 linkage groups consisting of 95 DNA 
markers. The average distance between markers was 7.0 
cM. Most of the linkage groups were consistent with 
those described previously [13]. 

QTL analyses were performed on the means of fatty 
acid composition, total oil content and expression levels 
of Olel, Olell and Olelll (Table 3; Figure 1). We 
detected 18 QTLs and 3 eQTLs for all traits examined. 
Individual eQTL or QTL were detected with percentage 
of variation explained (PVE or r^) 2.3% to 36.0%, and 5 
of them had PVE exceeding 10%. QTLs or eQTLs with 
positive and negative allelic effects were identified, with 
a positive effect implying a higher value for the trait 
conferred by the allele from PZM16 and vice versa (Fig- 
ure 2). 

QTLs for fatty acid composition and total oil content 

Eighteen QTLs were identified dispersed among all the 
linkage groups except LGs 3 and 11. A QTL of highly 
significant effect was determined to be located on LGl 
explaining 36% of variation of CI 8:1 composition, and 
was found to be associated with C18:2 compositions 
(Figure 2). Interestingly, another QTL on LGIO 
explained 5.9% of variation of C18:l composition was 
also associated with CI 8:2 compositions. Higher values 
for C18:l were conferred by the allele from PZM16, 
while higher values for CI 8:2 from Hybrid CI7041. 



Four QTLs were detected underlying total oil content. 
At the three QTLs on LGs 1, 2 and 4 respectively, the 
alleles from hybrid CI7041 contributed high total oil 
content. The most effective QTL was spotted on LG4 
explaining 11.1% of the variation, whose higher value for 
total oil content was conferred by the allele from hybrid 
CI7041. 

Favorite allele's effects 

There were strong QTLs for C18:l and total oil content 
detected on LGs 1 and 4, respectively. Mean phenotypic 
values of each trait were calculated for those progeny 
with the alternate alleles of the microsatellite markers, 
inherited from the /. integerrima (aa) or /. curcas (AA). 
A two-way ANOVA was performed on the progeny 
using two allelic combinations (AA, Aa) from markers 
linked to QTLs in order to investigate associations 
between phenotypic traits and genotypes of the QTLs. 
The phenotype values of each allelic combination of the 
QTLs are listed in Figure 3. Significant differences of 
phenotype means among different allelic combinations 
were identified, revealing the effects of alternative alleles 
inherited from the parents. 

Progenies with AA genotype at the marker Jcuint057 
located in qC18:l-l, showed the higher C18:l content 
(43.0%) than Aa (30.9%). By contrast, progeny with Aa 
genotype at the marker Jatr872 located in qOilC-4, 
showed the higher total oil content (38.0%) than AA 
(33.7%) (Figure 3). These results suggested the effect of 



Table 2 Correlation coefficients and significance of correlations among fatty acid composition, total oil content, 
oleosin gene expressions in a jatropha QTL mapping population 



Traits 


CI 6:0 


CI 8:0 


C18:1 


CI 8:2 


Total oil content 


Olel expression 


Olell expression 


CI 8:0 


-0.147 














C18:l 


0.038 


0.05 












CI 8:2 


-0.270** 


-0.155 


-0.962** 










Total oil content 


-0.087 


-0.02 


-0.003 


0.028 








Olel expression 


0.216** 


-0.029 


-0.013 


-0.038 


0.161* 






Olell expression 


0.139 


-0.047 


-0.074 


0.043 


0.191* 


0.696** 




Olelll expression 


0.132 


-0.005 


0.136 


-0.157 


0.170* 


0.790** 


0.697** 



P values are as follows: P < 0.05, P < 0.01. 



Liu et al. BMC Plant Biology 201 1 , 1 1 :1 32 
http://www.biomedcentral.eom/1471-2229/1 1/132 



Page 4 of 9 



Table 3 QTLs for seed oil traits and eQTLs for Olel, Olell and Olelll expressions in jatropha 



Trait 


QTL^ 


Linkage 


Marker 


Position'' 


LOD 




Additive 




(eQTL) 


Group 




cM 


Peak 


(%) 


Effects^ 


C16:0 (%) 


qC16:0-2 


2 


Jcuintl43 


474 


2.6 


0.1 


1.36 




qC16:0-7 


7 


Jatr802 


52.1 


3.1 


74 


1.42 




qC16:0-9 


9 


Jatr859 


15 


2.6 


7.2 


1.39 


C18:0 (%) 


qC18:0-2 


2 


curcin2 


52.6 


2.6 


5.3 


-0.69 




qC 18:0-5 


5 


Jatr746 


37.3 


6.9 


13 


1.15 




qC18:0-6 


6 


Jcuint036 


64 


3.9 


7.1 


-0.84 




qC 18:0-7 


7 


Jatr883 


40.3 


2.3 


4 


0.59 




qC18:0-9 


9 


Jatr859 


0 


9.2 


17.9 


1.26 


08:1 (%) 


qC18:l-l 


1 


Jcuint057 


0 


184 


36 


11.69 




qC18:l-5 


5 


Jatr739 


45.1 


2.3 


34 


-3.77 




qC18:l-10 


10 


Jcuintl80 


15.2 


4 


5.9 


4.75 


C18:2 (%) 


qCl 8:2-1 


1 


Jcuint057 


0 


16.5 


34.1 


-12.07 




qC18:2-6 


6 


Jatr301 


15 


24 


3.8 


4.3 




qCl 8:2-10 


10 


Jcuintl80 


15.2 


3 


4.6 


-4.4 


Total oil content (%) 


qOilC-1 


1 


Jatr722 


55.1 


2.3 


4.6 


-3.72 




qOilC-2 


2 


Jcuintl43 


474 


2.5 


4.9 


-3.74 




qOilC-4 


4 


Jatr872 


29.6 


5 


11.1 


-5.56 




qOilC-9 


9 


Jatr698 


18.6 


2.5 


5.2 


3.74 


Olel expression (AACt) 


qOlel-8 


8 


Jcuint277 


58.2 


1.9 


5.3 


1.71 


Olell expression (AAQ) 


qOlell-6 


6 


Jatrl52 


934 


2.6 


64 


-2.38 


Olelll expression (AAQ) 


qOlelll-5 


5 


Jatr739 


46.2 


3.1 


11.7 


-3.06 



^ QTL (eQTL): starting with "q," followed by an abbreviation of the trait name, the name of the linkage group, and the number of QTLs (eQTLs) affecting the trait 

on the linkage group. Olel, Olel expression level; Olell, Olell expression level; CI 6:0, CI 8:0,1 8:1 and CI 8:2, fatty acid compositions of palmitic acid (CI 6:0), stearic 

acid (CI 8:0), oleic acid (CI 8:1) and Linoleic acid (CI 8:2); OilC: Total oil content 

^ Position from the first marker on each linkage group. 

^ Proportion of phenotypic variance (R^) explained by a QTL (eQTL). 

^ Estimated phenotypic effect of substituting J. integerrima alleles with J. curcas alleles at QTL (eQTL). 



the two QTLs are opposite on these two key oil traits 
and favorite alleles were differentially from /. curcas and 
/. integerrima, 

eQTLs for oleosin genes 

SNP markers were developed in Olel, Olell and Olelll 
genes (Table 4), which were mapped on LGs 5, 3 and 5 
respectively (Figure 2). 

Olel and Olelll were mapped on LGS where the QTLs 
qCl8:0-5, qCl8:l-5 and qOlelll underlying C18:0, C18:l 
and Olelll expression clustered. Negative additive effect 
value of qOlelll-S indicated that /. curcas alleles were 
positive for Olelll expressions, of which LOD score was 
3.1. This eQTL of Olelll was localized near Olelll gene 
and overlapped with the QTLs controlling C18:0 and 
C18:l, revealing a cis- or trans-element for Olelll which 
subsequently controlling the C18:0 and C18:l. 

One eQTL on LG8 qOlel-8 was detected underlying 
Olel expression with LOD 1.9 (Table 3; Figures 1 and 
2). Additive effect value of qOlel-8 was positive, indicat- 
ing that /. integerrima alleles were positive for Olel 
expressions. To find as many putative QTLs (eQTLs) as 



possible, and to obtain a clearer understanding of the 
relationships among examined traits, a threshold eQTL 
of 1.9 for declaring a suggestive eQTL was employed. 
Low thresholds may not be useful in plant breeding pro- 
grams but they have been shown to help in understand- 
ing relationships among traits [18]. 

Olell was located on LG3. One eQTL for Olell was 
detected on LG6 with LOD 2.6, which closed to qC18:0- 
6. It is suggested that a trans-element for Olell could 
harbor in this region which controlling the C18:0. Addi- 
tive effect values indicated that /. curcas alleles were 
negative, indicating that the effect of /. curcas alleles was 
positive for Olell expressions. 

Discussion 

Development of inter-specific populations 

To broaden the genetic diversity of cultivated crops and 
to identify QTLs associated with beneficial traits, such 
as yield, grain quality and disease resistance, develop- 
ment of inter-specific populations is a feasible strategy 
[19]. We developed around 500 SSR markers in jatro- 
pha, but very low polymorphism was detected within /. 
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Figure 1 Whole genome scan for QTL for oil traits and Oleosin 
gene expression in jatropha. A QTL scans of oil traits on linkage 
maps. Horizontal line indicates 5% LOD significance thresholds (2.5) 
based on permutation. B QTL scans of Olel, Olell and Olelll 
expressions on linkage maps. Horizontal line indicates LOD 
significance threshold (2.0). 



curcas, indicating the genetic variation was very limited 
within /. curcas. Thereby, we successful constructed a 
QTL/eQTL mapping population by crossing /. curcas to 
/. integerrima and generating a backcrossing population, 
and observed enhanced genetic diversity on DNA, RNA 
and phenotype levels, which was the prerequisite for 
QTL and eQTL detection. 

For oil trait improvement, the interspecific hybridiza- 
tion approach is also viewed as a viable method to intro- 
gress the traits of interest, i.e. namely more liquid olein 
in oil palm [20]. With MAS, selection can be carried 
out in segregating generations of interspecific hybrids 
and their backcrosses more discriminately using molecu- 
lar markers linked to the specific fatty acids. We investi- 
gated effects of the QTLs on oil traits and found that 
favorite alleles were originated from not only /. curcas 
but also /. integerrima, CI 8:1 in /. curcas was higher 
than in /. integerrima, while total oil content in /. inte- 
gerrima was 51.04%, much higher than in /. curcas 
(Table 1). Consistent to this result, qC18:l-l and qOilC- 
4, controlling C18:l and total oil content respectively, 
were detected with the favorite alleles originated from /. 



curcas and /. integerrima respectively. Therefore, the 
QTL mapping population will be very useful for trans- 
ferring favorite alleles form the two parents by further 
backcrossing and marker assisted selection. 

Various germplasms were successfully utilized for 
development of chromosome segment substitution lines 
for studies on pest and disease resistance and other 
agronomic triats in rice [21-23]. Here we generated 
backcross populations for map construction and QTL 
mapping, which required less time to be developed and 
being 'immortal' for future QTL mapping due to jatro- 
pha's perennial life cycle. Meanwhile, the specific advan- 
tage of backcross populations is that, the populations 
can be further utilized to develop chromosome segment 
substitution lines for marker-assisted backcross breed- 
ing. The chromosome segment substitution lines will 
provide a valuable tool for jatropha germplasm enhance- 
ment, and can be expected to reveal the genetic basis of 
traits specific to the donor /. integerrima. 

Linkage or pleiotropic effect of genes in QTL cluster 

Several chromosomal regions were associated with more 
than two traits indicating either linkage or pleiotropic 
effect. We detected a QTL cluster controlling CI 8:1 and 
CI 8:2 contents on the same region, i.e. closed to marker 
JcuintOS? on LGl and Jcuintl80 on LGIO with the addi- 
tive value of C18:l opposite to that of C18:2. This could 
explain the strong negative correlation between CI 8:1 
and C18:2 (Table 2), which was consistent to the fact 
that linoleic acid is desaturated from oleic acid. Espe- 
cially on LGl, the QTL was detected with a highly sig- 
nificant effect, accounting for 36.0% of the variation. It 
is revealed that either certain genes coexisted in these 
QTLs or a certain gene with pleiotropic effect in fatty 
acid metabolism pathway by modulating both CI 8:1 and 
C18:2 contents simultaneously. It will be meaningful to 
conduct fine mapping of these QTLs, isolate the target 
genes, and understand whether linkage or pleiotropic 
effect. The QTL regions were still distant to the flanking 
markers with linkage distance larger than 2 cM. Fine 
mapped QTL will speed up genetic improvement 
through MAS [3]. Construction of a high-resolution 
genetic linkage map of jatropha is underway, which will 
lay a solid foundation for a variety of future genetic and 
genomic studies, including QTL fine mapping and mar- 
ker assisted selection. 

eQTL analysis of oleosin genes 

To examine the function and modulation of oleosin 
genes in jatropha, we determined the expression levels 
of Olel, Olell and Olelll in the QTL mapping popula- 
tion, and conducted analysis with an approach named 
"genetical genomics" for identifying the genomic regions 
influencing gene expression [12,24]. The correlation of a 
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Figure 2 Summary of QTL (eQTL) locations detected on thie genome of jatrophia. QTLs (eQTLs) represented by bars are shown on the left 
of the linkage groups, close to their corresponding markers. The lengths of the bars are proportional to the confidence intervals of the 
corresponding QTLs (eQTLs) in which the inner line indicates position of maximum LOD score. 
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structural gene's map position and its eQTL provides an 
indication of its regulation [24]. If the position of one 
gene and its eQTL are congruent, cis-regulation could 
be inferred, which means that the allelic polymorphism 
of the gene itself, or closely linked regulatory elements, 
directly impact the gene's expression. In this study, the 
eQTL for oleosin genes do not colocalize with these 
gene. This result suggests that the observed differences 
in oleosin gene expressions could be the consequences 
of trans- regulation, which means that gene expression is 
mainly regulated by trans-acting factors. A similar phe- 
nomenon has been observed for a set of genes involved 
in the biosynthesis of lignin in Eucalyptus. Most of 
these genes were significantly influenced by two eQTLs 
on LGs 4 and 9, whereas the structural genes were dis- 
tributed throughout the entire genome [25]. 

The significant but low correlations were observed 
between oil traits and the expressions of oleosin genes. 
Similar phenomenon was reported by Yin et al [10]. 
They reported that the significant correlation between 
the expression of both GmRCAa and GmRCAb and 
Rubisco initial activity, photosynthetic rate, and seed 
yield indicated that these genes could play a role in 
increasing photosynthetic capacity and seed yield. How- 
ever, the correlation coefficients between gene expres- 
sion and Rubisco initial activity, photosynthetic rate, and 
seed yield were relatively small. This was also reflected 
by the fact that no coincident QTL (eQTL) was found 
between gene expression levels and the other three 
traits. Thus, they concluded that factors other than 
GmRCAa and GmRCAb limited photosynthetic capacity 
and seed yield. In our study, significant but low correla- 
tions between oil traits and the expressions of oleosin 
genes indicated that these genes could affect fatty acid 
composition and content; meanwhile, there should be 
other complex factors together with oleosin genes affect- 
ing oil traits. 

The three eQTLs will provide possible approaches to 
oil trait improvement beyond previous QTL mapping 
results. Interestingly, Olelll gene, eQTL of Olelll 
qOleIII-5 and QTL of qC18:0-5 and qC18:l-5 were clus- 
tered on the same region on LGS. To further address 
whether a cis- or trans-element for Olelll harbored on 
LGS subsequently controls the fatty acid compositions, 
fine mapping the two loci is still needed. Only eQTL of 
Olelll was coincident with QTL for oil composition, this 
result could be resulted from function differentiation of 
oleosin genes. 

Conclusions 

In conclusion, we identified 18 QTLs underlying the oil 
traits and 3 eQTLs of the oleosin acid genes. Among 
them, qC18:l-l, qOilC-4 and qOlelll-S, controlling oleic 
acid, total oil content and oleosin gene expression 



respectively, were detected with relatively high contribu- 
tion rates (R^) and could be expected to be applied in 
MAS by integrating more markers in these region. 
These data represents the first successful detection of 
QTLs/eQTLs underlying key agronomic traits in 
jatropha. 

Methods 

Plant material and plant growth conditions 

/. curcas PZM16 was crossed to /. integerima SOOl and 
hybrids CI7041 were generated. Then a backcrossing 
(BC) population was constructed consisting 286 indivi- 
duals derived from the backcross PZM16 x CI7041. The 
population and parental lines were planted under stan- 
dard growth conditions in experimental field of Lim 
Chu Kang farm, Singapore. 

Isolation of genomic DNA and synthesis of cDNA 

Total DNA from leaves was extracted and purified using 
the DNeasy plant mini kit (QIAGEN, Germany). Oil 
bodies are located inside the cells of mature seeds. Total 
oil content and fatty acid composition in mature seeds 
are agronomic traits of importance. To investigate 
expressions of oleosin genes in mature seeds which are 
used for oil extraction, total RNA was isolated from 
mature seeds using plant RNA purification reagent 
(Invitrogen). Poly(A) tails were then added to the 3' end 
of the RNAs by poly (A) polymerase (Ambion), and the 
polyadenylated RNAs were reverse transcribed by Super- 
Script II reverse transcriptase (Invitrogen) with the oligo 
(dT) 3'-RACE adaptor (Ambion). 

Trait measurement and data collection 

Each sample of QTL mapping population was grinded 
with liquid nitrogen, divided into 3 copies. Every sample 
consists of 3 mature seeds collected randomly from the 
same tree. Fatty acid compositions were analyzed by 
Gas chromatography (GC). Total lipid, extracted from 
100 mg mature seeds, was transmethylated with 3 N 
methanolic-HCl (Sigma, St. Louis, MO, USA) plus 400 
(iL 2,2,-dimethoxypropane (Sigma, St. Louis, MO, USA). 
Oil was extracted using solvent (hexane) extraction fol- 
lowed by esterification to transfer from oil to methyl 
ester. The fatty acid methyl esters (FAME) was analyzed 
by GC using GC Agilent 6890 (Palo Alto, CA, USA) 
employing helium as the carrier gas and DB-23 columns 
for components separation. The GC analytical method 
was performed at 140°C for SO s and a 30°C min"^ ramp 
to 240°C, and the final temperature was maintained for 
SO s for a total run time of 32 min. FA composition 
value included in the analyses was calculated based on 
peak area. 

To amplify the mRNA from the reverse transcribed 
cDNAs and determine expression levels, real-time PCR 
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was conducted with Real-Time PGR machine (I-Cycle, 
BioRad). Each reaction contained 200 ng of first-strand 
cDNAs, 0.5 [iL of 10 mmol L'^ gene-specific primers, and 
12.5 [iL of real-time PGR SYBR MIX (iQ™ SYBR® Green 
Supermix, Bio-Rad). Amplification conditions were 95°G 
for 5 min followed by 40 cycles of 95°G for 15 s and 60°G 
for 60 s. The jatropha 18S rRNA was selected as the 
endogenous reference was used as a control to test for 
sample-to-sample variation in the amount of cDNA. 
cDNA from mature seeds of jatropha PZM16 was used 
as the calibrator on each real-time PGR plate. Two tech- 
nical replicates of each reaction were performed. Nor- 
malized expression for each line was calculated as 

described in [10], i.e. AAGt = (Gt, Target - Ct, iss) genotype 

- (Ct, Target - Ct, i8s)caiibrator. Lower AAGt value means 
stronger gene expression and vice versa. Five mature 
seeds from each plant of QTL mapping population were 
used to determine the relative expression levels of Olel, 
Olell and OlellL The results presented are means of the 
biological replicates for each plant. 

DNA markers and genotyping 

Ninety- five markers almost evenly covering the 11 LGs 
were selected from a first-generation linkage map of 
jatropha [13]. One primer of each pair was labeled 
with FAM or HEX fluorescent dyes at the 5'end. The 
PGR program for microsatellite amplifications on PTG- 
100 PGR machines (MJ Research, GA, USA) consisted 
of the following steps: 94°G for 2 min followed by 37 
cycles of 94°G for 30 s, 55°G for 30 s and 72°G for 45 
s, then a final step of 72°G for 5 min. Each PGR reac- 
tion consisted of Ix PGR buffer (Finnzymes, Espoo, 
Finland) with 1.5 mM MgG^, 200 nM of each PGR 
primer, 50 (iM of each dNTP, 10 ng genomic DNA 
and one unit of DNA-polymerase (Finnzymes, Espoo, 
Finland). Products were analyzed using a DNA sequen- 
cer ABI3730xl, and fragment sizes were determined 
against the size standard ROX-500 (Applied Biosys- 
tems, GA, USA) with software GeneMapper V3.5 
(Applied Biosystems, GA, USA) as described previously 
[26]. 



Statistical analysis and QTL (eQTL) mapping 

QTL (eQTL) analysis allows the genetic basis of variation 
of quantitative traits of interest to be dissected. Scoring 
every individual of a mapping population for the trait of 
interest and establishing a genetic linkage map for that 
population are two prerequisites for QTL (eQTL) detec- 
tion. In this study, expression level data of fatty acid com- 
position and content, and Olel, Olell and Olelll expression 
levels of the backcross population consisting of 286 indivi- 
duals were collected with 3 replications. Pearson phenoty- 
pic correlations among traits were calculated by SAS 
PROG GORR. The 95 markers were genotyped in the 
QTL mapping population. SNP markers for mapping the 
three genes and primer pairs for determining expression 
levels by real-time PGR were listed in Table 4. 

Linkage map was constructed using the software 
GRIMAP 3.0 to detect Unkage and build map [27]. All 
multipoint distances were calculated using the Kosambi 
function. MapGhart 2.2 software was used for graphical 
visualization of the linkage groups [28]. QTL (eQTL) 
analysis was performed using QTL Gartographer version 
2.5 [29]. Model 6 of composite interval mapping was 
deployed for mapping QTLs (eQTLs) and estimating 
their effects. The genome was scanned at 2 cM intervals, 
and the forward regression method was selected. The log 
of the odds (LOD) score for declaring a significant QTL 
(eQTL) by permutation test analyses (1,000 permuta- 
tions, 5% overall error level) as described previously. To 
find as many putative QTLs (eQTLs) as possible, and to 
obtain a clearer understanding of the relationships 
among examined traits, a threshold eQTL analysis of 
oleosin genes in of 2.0 for declaring a QTL (eQTL) was 
employed. Low thresholds may not be useful in plant 
breeding programs but they have been shown to help in 
understanding relationships among traits [18]. 

The maximum LOD score along the interval was 
taken as the position of the QTL (eQTL), and the region 
in the LOD score within 1 LOD unit of maximum was 
taken as the confidence interval. Additive effects of QTL 
(eQTL) detected were estimated from composite interval 
mapping results as the mean effect of replacing hybrid 



Table 4 SNP markers and real time PCR primer pairs for Olel, Olell and Olelll genes 


Gene 


Forward primer (5'-3') 
Reverse primer (5'-3') 


PCR product length 
(bp) 


For SNP or Real time PCR use 


Olel 


CATOCGCTAGCTGTOCGACTCC 
CGCCGCmGCCAmCCATCT 


207 


SNP and Real time PCR 


Olell 


GGGGCTATGGGGCTCACAG 
GTOAGTOGmATGGGGGATCT 


313 


SNP and Real time PCR 


Olelll 


ACAGCCACGATCCCACCAAGTAGT 
GGACAGAGCTGAGCAGmGGACA 


443 


SNP 


Olelll 


tggtgccgacgg™tcac 
tacatgctgtccaaactgctcag 


216 


Real time PCR 
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(CI7041)'s alleles at the locus of interest by /. curcas 
(PZM16) alleles. Thus, at a QTL (eQTL) having a posi- 
tive effect, the alleles of /. curcas will increase the trait 
value. The contribution of each identified QTL (eQTL) 
to total phenotypic variance (r^) was estimated by var- 
iance component analysis. QTL (eQTL) nomenclature 
was adapted as following: starting with "q," followed by 
an abbreviation of the trait name, the name of the link- 
age group and the number of QTL (eQTL) affecting the 
trait on the linkage group. 

In order to investigate associations between phenotypic 
traits and genotypes of two QTLs on LGs 1 and 4, mean 
phenotypic values of traits were calculated for those pro- 
geny with the alternate alleles of the microsatellite mar- 
kers, inherited from the /. integerrima (aa), alleles 
inherited from the /. curcas (AA). A two-way ANOVA 
was performed on the progeny using two allelic combina- 
tions (AA, Aa) from markers linked to QTLs. This was 
conducted by using the general linear model (GLM) pro- 
cedure of SAS (SAS Institute) and the Bonferroni method 
of multiple comparisons with a < 0.01. 
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